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Abstract 


In multidisciplinary optimization problems, response surface techniques can be used to replace the 
complex analyses that define the objective function and/or constraints with simple functions, typically 
polynomials. In this work a response surface is applied to the design optimization of a helicopter rotor 
blade. In previous work, this problem has been formulated with a multilevel approach. Here, the 
response surface takes advantage of this decomposition and is used to replace the lower level, a 
structural optimization of the blade. Problems that were encountered and important considerations in 
applying the response surface are discussed. Preliminary results are also presented that illustrate the 

benefits of using the response surface. 


Introduction 

In dealing with multidisciplinary design optimization (MDO), several problems arise. The analysis 
codes of the various disciplines are vastly different and difficult to integrate. In addition, these analyses 
must be linked to an optimizer, which is another difficult task. Because of their complexity, 
multidisciplinary problems are typically decomposed into several simpler subproblems. This process is 
known as multilevel decomposition (Refs. 1-6), and the functions that pass information between the 
various levels may not always be smooth. Moreover, individual disciplines may require more than one 
computational model. For example, the finite-element structural analysis may employ a simplified 
aerodynamic code to predict the loads, and the Computational Fluid Dynamics (CFD) analysis may use 
a simplified structural model such as the weight equation. Researchers are investigating approximation 
concepts to reduce the high computational costs associated with MDO problems. Barthelemy ( Ref > 
reviews some of the most recent approximation concepts used in structural optimization problems. As 
discussed in Ref. 7, one area which is gaining attention is response surface methodology (e.g., Rets. / 

15). 
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A response surface is a means of approximating complex functions and/or analyses with a 
polynomial function or a neural network. Response surfaces will be discussed in more detail later in this 
paper. The benefits of a response surface approach to MDO are three-fold. First, the response surface 
provides a method for smoothing the noisy functions (Refs. 7 and 8). Eason and Fenton (Ref. 9) suggest 
that many complex engineering analyses (e.g., finite-element models) yield nonsmooth functions that 
may cause errors in the finite difference derivative calculations commonly used in conventional 
optimization algorithms. Thus, these analyses should be replaced with smooth analytical functions. 
Second, the response surface approximation provides one discipline with a quick assessment of other 
disciplinary analyses, similar to empirical structural weight equations. Third, response surfaces are 
amenable to parallel computation, whereby the various disciplinary analyses can be conducted on 
different processors or even different machines (Ref. 7). 

References 7-14 describe applications where response surfaces have been used in structural 
optimization problems. In these applications, response surfaces are generated to develop analytical 
functions in place of the complex analyses or experimental data. These analytical functions are then 
used in the optimization procedure. The work presented in this paper applies response surface 
methodology to a rotor blade optimization procedure (Ref. 16) which used multilevel decomposition. 
Unlike the applications in Refs. 8-14, a response surface is generated to approximate an entire 
optimization procedure. Specifically, the lower level structural optimization in Ref. 16 is replaced by a 
response surface. First, a summary of the multilevel rotor blade optimization procedure (Ref. 16) is 
given. Next, a discussion of response surface implementation is discussed. Finally, results are 

presented. 


Integrated Aerodynamic/Dynamic/Structural (IADS) Optimization Procedure 

Reference 16 describes an integrated aerodynamic/dynamic/structural (IADS) optimization 
procedure for helicopter rotor blades. The procedure combines performance, dynamics, and structural 
analyses with a general purpose optimizer (CONMIN, Ref. 17) using multilevel decomposition 
techniques. The upper level optimizes the blade by changing global quantities such as blade planform, 
twist, and stiffness distributions. The upper level chord and stiffness distributions are treated as 
independent quantities. The reconciliation between these distributions is done on the lower level, which 
consists of several independent subproblems at stations along the blade radius. These subproblems 
optimize detailed cross-sectional dimensions to assure structural integrity of the blade and to reconcile 
the upper level independent chord and stiffness distributions with the lower level calculated stiffness 
distributions. This reconciliation is improved further by a set of upper level coordination constraints. 

Upper Level Optimization 

For the present work, the blade is assumed to have a rectangular planform and constant stiffnesses 
spanwise. Thus, a subset of the upper level design variables described in Ref. 16 is used. These design 
variables (Fig. 1) are the maximum pretwist 0 tw , the root chord c r and the four blade stiffnesses: 
chordwise bending E^, flapwise bending EI^, extensional EA, and torsional GJ stiffness. 
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Figure 1. Upper level design variables. 


The upper level objective function to be minimized is a measure of the horsepowers required in each of 
teX* conditions (hover, forward flight, and maneuver) and the hub shear (t.e„ vtbratron) and ts 

given by 
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where P h P ff and P m are the powers required in hover, forward flight, and maneuver, respectively. N is 
th^number of bladesTand S Nff is the N per rev rotating vertical hub shear in forward flight. The terms 

ki , ko, k 3 , and k4 are weighting factors chosen by the user. Ph ref * Pff„f • ^ m ref ’ anc ^ 316 
values used to normalize and nondimensionalize the objective function components. The constraints on 
tta prebTem arc aerodynamtc and dynamic in nature (see Ref. 16 for details) The coordmauon 
constraint which is the means of communication between the upper and lower levels is also tmposec ' a > 
the upper level The Langley-developed hover analysis HOVT (a strip-theory momentum analysis based 
used to predict the power required in hover. The comprehenstve rotor analysts 
CAMRAD/JA (Ref. 19) is used to predict forward flight and maneuver performance. 


Lower Level Optimization 

The purpose of the lower level is to assess 
to provide the stiffnesses required by the upper 


whether a structure at a given radial location can be sized 
level optimization and still have the strength to withstand 
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loads calculated by the upper level analysis. In the present work a single cross section is optimized. 
The lower level represents the structural properties of the blade in terms of wall thickness 
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Figure 2. Lower level design variables. 


and lumped areas in the boa beam cross section (Fig. 2). The lower level objective function to 
minimized is given by 
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and is a measure of the stiffness matching. The stiffness matching is the difference between the 
stiffnesses that are desired by the upper level (denoted by the superscript U) and the stiffnesses that can 
be obtained by the lower level (denoted by the superscript L). Constraints at the lower level assure 
structural integrity of the blade (see Ref. 16). For convenience, the set of lower level constraints is 
replaced by a single cumulative constraint, an envelope function known as the Kreisselmeir-Steinhauser 
(KS function. Ref. 20) function which approximates the active constraint boundary 


KS = g max +-ln| 


N 


X* 

j= 1 


P(gc, gmu ) 


<0 


(3) 


where g max is the maximum constraint component, nc is the number of lower level constraint 

components and p is defined by the user. The lower level optimization is repeated for various values of 
p. Initially p is small and then increases until a maximum value p,^ is reached. For large values ot p, 
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the value of KS approaches g max . The KS function is a single measure of the degree of constraint 
satisfaction or violation and is positive (violated) if at least one of the constraints §cj is violated. 

Coordination 

The coordination between upper and lower levels is implemented by an upper level constraint 
imposed to encourage changes in the upper level design variables which promote consistency between 
the upper and lower level stiffnesses. Specifically, this constraint has the form 

g = F U -(l + e)F o L pt <0 (4) 

where F^ pt i s the most recent value of the lower level objective function (i.e., optimum value of 

equation 2), F u is an estimate of the change in F opt caused by a change in the upper level design 
variable values, and e is a specified tolerance denoting the coordination parameter. The coordination 
constraint relays to the upper level the value of the lower level optimal objective function (i.e., how well 
the stiffnesses match) and how that value will change when an upper level design variable changes. 

Figure 3 shows a flowchart of the optimization procedure. In this research, the entire lower level 
structural optimization is replaced by a response surface. The surface is generated outside the 
optimization cycle loop and is used both in place of the lower level during the optimization process and 
in the coordination constraint. The implementation of the response surface in this problem will be 
discussed later. 



Figure 3. Flowchart of IADS optimization with response surface formulation. 
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Response Surface 


As previously mentioned, the response surface is an approximation that replaces the objective 
function and/or constraints with a polynomial expression or a neural network. The reader is referred to 
Refs. 7 and 15 for more details on response surface methodology. A response surface was constructed 
to replace the entire lower level structural optimization procedure in the form of the optimum lower 

level objective function (Fo pt ). Therefore, the response surface predicts the stiffness mismatch between 
the upper and lower levels. A quadratic polynomial of the form 

n n n 

Fopt =b 0 +Xb i x i +XX b ij x i x j (5) 

i=l i=l j=i 

was selected as the model. Here, the coefficients b are computed by a least-squares regression analysis. 
The response surface is a function of the upper level design variables \ v In this case, 6 design variables 
yield 28 undetermined coefficients for which to solve. The points in the design space used to generate 
the response surface are given by statistical experimental design (specifically, a central composite 
design). Figure 4 shows a central composite design in three dimensions. A set of baseline design 
variables (i.e., an initial guess for the optimizer) is scaled to the point (0, 0, 0, 0, 0, 0). The baseline 
design variables are perturbed by ±10 percent. Several larger percentages (e.g., ±25) were tried. The 
CAMRAD/JA analysis experienced numerical difficulties in analyzing some of the design points 
generated with larger perturbation values. For the work presented here, the ±10 percent was found to be 
the best value. The design variables are then scaled to ±1. These points form the vertices of the 
hypercube. Finally, the centers of each face of the hypercube are selected. The 6 design variables in 
this problem used 77 design points, which creates an overdetermined system. 



( 0 , 1 , 0 ) 


Figure 4. Central composite design in three dimensions. 


Important Considerations 


Response Surface Implementation 

In conducting this research, several important considerations must be addressed in the application 
of a response surface to this type of problem, in which an entire optimization is replaced with a 
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polynomial function. First one must determine what the response surface will approximate; two 
possible choices are the optimal lower level objective function value ( F opt) or the same objective 
penalized by the lower level constraints ( F o L p. + P®** KS function). A third possibility that was not 
explored is to build a second response surface for the lower level KS function CMdriji. ^and use rfhts 
approximation as an additional upper level constraint. Introduction of the constraints from the lower 
level ensures that the designs are feasible; however, for the results presente ere east y 
concern because the lower level constraints of the generated designs were sattsfted. Therefore, the 
unpenalized optimal lower level objective function was selected to be approximated by the response 

surface. 

CM “n"m.ion tha, requires attention is the selection of the coordinatton constraint. 
Recall that the coordination constraint is the means by which the upper and lower levels 
In essence, the coordination constraint discourages the upper level opt.nuzat.on from degrading .the 
stiffness matching achieved on the lower level. As formulated in Ref. 16, the lower level provides 

upper level not only with F o L P t, but also with the optimal values of the lower level design. F tn eqn. 4, 
is evaluated at the upper level using the current upper level design variable values and the previous 
optimum lower level design variables. F^ t is the previous optimum lower level objective function. 
When a response surface formulation is used, however, all dependence on lower level des.gn vanabtes is 
lost; the values of the lower level design variables are not available. The approxtma ing P° y 
only a function of the upper level design variables. One possible solution to tlus problem is i to bn 
response surface for each of the six lower level design variables, as well as for the lower level 
optimization. Another solution is to reformulate the coordination constraint as follows 

g = F<5i>t^Y(Y = 0.001) (6) 

This constraint forces the matching of stiffnesses to be within a small tolerance (?) of zero. Now the 
coordination constraint can be evaluated by the response surface (eqn. 5). 

Oualitv of Design Points Used to Generate the Response Surface 

The third issue that must be examined is whether or not the response surface represents an optimum 

design In the present work, the response surface is formulated as a quadratic approximation of F o P t. 
The coefficients of this polynomial are determined from the 77 structural optimizations used to generate 
the response surface. Therefore, the response surface approximation is only as good as those 77 design 
poin? If those designs are not optimal, then the response surface will be ineffective m predicting the 
stiffness mismatch (its desired task). Below is a situation encountered in this study where optimality of 

individual points was a problem. 

Initially, the response surface was generated from a set of 77 optimum lower level designs^ Each 
design was obtained using the same minimum and maximum values for p in the KS function (Eqn. ) 
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used in Ref. 16. When the structural optimization at the lower level was replaced with t is response 
surface, the upper level (overall) optimization consistently converged to the set of baseline design 
variables given as an initial guess, which was not an optimum. To assess the problem, the generation o 
the response surface data points was examined. The lower level structural optimization was performed 
with different initial values of wall thicknesses. (Note: initially the lumped areas are set to zero.) Each 
time the initial thickness values were changed, the lower level objective function converged to a 

different value These different values suggest the possibility of many local optima in the lower 
level objective function. However, this possibility was discovered not to be the case; a test was 
conducted in which the objective function was evaluated at regular intervals between two optimum 
points (a, and x 2 ). These intermediate points were governed by a parameter a (0 < a < 1) in the 

following expression: 


x = ax, + (l-a)x. 


m 


Several pairs of points were tested; a sample of the results is presented in Fig. 5. This plot shows that 
the lower level objective function does not experience any local maxima or minima. 



Figure 5. Test for local optima in lower level objective function. 

Next the maximum value for p was increased. Then the structural optimizations at the lower level 
for all initial thicknesses yielded the same lower level optimum. The response surface was then 
regenerated with 77 new structural optimizations and the overall optimization no longer converged to 
the set of baseline design variables. The reason that a global optimum was not found was attributed to 
the value of p in the KS function (Eq. (3)). The maximum p value was too small; thus, the structural 
optimization converged prematurely. This response surface was used to compute the results presented 

in this paper. 
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'^Sy’o'L the response surface has been generated, the qual.ty of the fit must be assessed Two 
approaches were considered for this assessment. The first approach is to measure the percent error 
between the response surface and the actual structural optimization 


% error = 100 


pL _ pL 
r opt r opt 


A n f thp prrnr is nossible at individual points in the design space. However, as discussed 

earfc Zty mtJevel functions can be noisy which is illustnned by Ftg. 6. Thus, although individua^ 
pohns may not be Cose to the response surface values, the overall fit of £ 
adequate. Carpenter (Ref. 14) suggests a nondimenstonal root mean square (RMS) error as 

this overall fit 

t(F 0 L pt, -FopO 2 

— n — 

RMS error = 100 ^ ^ 


Here, F^ t is the average 

F^t is s iven b y 


value of the optimal lower level objective function over N points examined; 
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Figure 6. Response surface approximation. 
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The RMS error computed for the set of 77 points used to generate the response surface was 1 8 
percent. The least-squares regression analysis that computes the coefficients of the approximating 
polynomial seeks to minimize error at these points. A better test of the quality of fit, therefore, is to 
compute the RMS error at 50 random points in the design space; this value was 59 percent. Both errors 
are reasonable and are comparable to the results reported by Carpenter (Ref. 14). A total of 127 pomts 
were used to generated and validate the response surface. 


Results 

Both the response surface formulation and the IADS formulation with the modified coordination 
constraint (Eq. (6)) were studied. Both optimization routines were implemented on the NASA Langley 
Cray 2 supercomputer. Results are presented using the parameters in Table 1. Table 2 presents the 
upper level design variables and objective function for the initial blade design. The upper level 
objective function OBJ is also given in the table. This design is infeasible because an upper level 
frequency constraint is violated. Neither formulation was run to convergence. Results are presented 
after 15 cycles and 30 cycles. The response surface was only generated once. 


Table 1. Parameters 


N, number of blades 

k l 

k 2 

k 3 

k4 

Ph ref 

Pff tef 

Pm ref 

S4 ref 

4 

10.0 
5.0 
5.0 
0.5 
15 hp 

13 hp 
12 hp 
2 lbf 

Table 2. Initial Guess for Blade Design 

Upper Level Design Variable 

Initial Design 

Twist (deg) 

-9.00 

Chord (ft) 

0.4500 

EI XX (lb-ft 2 ) 

2,907 

EI ZZ (lb-ft 2 ) 

226.6 

GJ (lb-ft 2 ) 

261.9 

EA (lb) 

1,466,000 

OBJ 

20.89 
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Table 3 compares the results of the two formulations after 15 optimization cycles. In both cases, 
the value of OBJ has decreased from the initial value. The design variable values are similar, except for 
the chord value. At first glance, the IADS formulation seems to perform better than the response 
surface. However, one also needs to examine what is happening at the lower level (i.e„ the sttffne 
matching). Recall with the response surface formulation, the actual wall thicknesses and lumpe areas 
are not available. To determine these dimensions and to assure the structural integrity of the response 
surface generated blade design, an lower level structural optimization was performed at the lower level 
to determine how well the stiffnesses match. 


T^l^^^rimizationResult^fte^l5 i C^clM 


Upper Level Design Variable Response Surface 

i ^ i rv 


Twist (deg) 
Chord (ft) 

EIxx Ob-ft 2 ) 

EIzz (lb-ft 2 ) 

GJ (lb-ft 2 ) 
EA (lb) 


- 12.10 
0.4381 
2,831 

231.9 

225.4 

1,469,000 


IADS 

-12.55 

0.3553 

2,731 

212.5 

225.0 

1,471,000 

19.98 


OBJ 

p L 

r opt 


20.87 

0.004 


0.5474 


The upper level design variables were taken at the end of the 15th cycle (using the response surface) and 
dTset of thicknesses and areas that gives these upper level stiffnesses was found. No approximation 
was made here. Figure 7 shows the stiffness matching between the upper and lower levels, the response 
surface (shaded) is compared with the IADS formulation (black). The response surface design after 15 

cycles matches the stiffness well ( Fjj* = 0.004). The IADS design matches the stiffnesses poorly ( opt 
= 0 5474) These plots indicate that, although OBJ for the IADS procedure is smaller than that found by 
the response surface procedure, the design is infeasible because of poor matching, l«todarly 'in .the 
flapwise bending and torsional stiffnesses. The response surface design, nevertheless, still violates 
coordination constraint, which states that the matching must be within 0.001 (Ti- 
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Figure 7. Stiffness Matching After 15 Cycles 


For the next 15 cycles, the move limits in the IADS formulation were cut in half in order to improve 
the results. Table 4 shows the optimization results after 30 cycles. Stiffness matching at the end of the 
30th cycle is shown in Fig. 8. As seen in the plots, the matching of the IADS procedure has improved 
only slightly; the result remains an infeasible design. The smaller move limits, therefore, had little 
effect on the performance of the IADS procedure. On the other hand, for the response surface the 

matching has improved and F 0 L pt has decreased to 0.0005; thus this design is feasible. 
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Twist (deg) 
Chord (ft) 
EIxx Ob-ft 2 ) 
Elzz (lb-ft 2 ) 

GJ (lb-ft 2 ) 
EA (lb) 

OBJ 


-10.55 
0.4328 
2,826 

210.5 

228.0 

1,470,000 

20.73 


- 1 1 .< 
0.3504 
2,800 

214.3 

228.0 

1,525,000 


0.0005 


19.93 

0.6328 


“opt 


Response surface formulation 



IADS formulation 
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Figure 8. Stiffness matching after 30 cycles. 
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_ . n __ : n thp move limits for the IADS formulation resulted in no apparent improvement 

response surface resul^ The has improved significantly, but 

» that the design ,s in^He. H o.eve, te e ■ A DS 
results are similar to those obtained with the response surface. 


1 Response surface formulation 

■■ IADS formulation (with originial coordination constraint) 
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Figure 9. Stiffness matching After 15 cycles 
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Table 5. Optimization Results After 30 Cycles 
(IADS formulatioirusesorigndco ^in^on^con^^m^^n^ 
Upper Level Design Variables Response Surface IADS 


Twist (deg) 
Chord (ft) 
EI XX (lb-ft 2 ) 

EIzz (lb-ft 2 ) 

GJ (lb-ft 2 ) 
EA (lb) 

OBJ 


- 12.10 
0.4381 
2,831 

231.9 
225.4 
1,469,000 
20.67 


-11.49 

0.3836 

3,011 

163.0 

253.3 

1,482,000 

20.23 


NtI . , he resDOnse formulation and the IADS formulation are examined on the basis of the number 
Next the response response surface, 127 structural optimizations were 

of structural optimizations requned. In b “ ldm 8 * P sit( . design plus the 50 random points for 

formulation to determine the structur cvc les than the IADS formulation requires, 

formulation requires more structural optimizations aft y , i desien The 

However, the response formulation design is much ^Scyclel, 

opJSizations and the response 

” * 4 .* tu T a nc formulation still has poor stiffness matching. The 

additional structural optimizations. The IADS formulation su v typically 

:rs - - - • - • 

design. 


Concluding Remarks 

f -- hoc hf*pn developed to replace the lower level structural optimization of a 
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panted in this paper are encouraging and confirm that response surfaces can be a vtable way of eas g 
the computational demands of multidisciplinary design optimization. 


MgM • The research of the first Mho. was luridly supported utrier the Ungley Aerospace Research Summer 
Scholars(LARSS) Program at NASA Langley Research Center in Hampton. Virginia. 
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